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ABSTRACT 

We present observations of the radio emission from the unusual supernova SN 1988Z in MCG +03- 
28-022 made with the Very Large Array at 20, 6, 3.6, and 2 cm, including new observations from 1989 
December 21, 385 days after the optically estimated explosion date, through 2001 January 25, 4,438 days 
after explosion. At a redshift z — 0.022 for the parent galaxy (~100 Mpc for H = 65 km s _1 Mpc -1 ), 
SN 1988Z is the most distant radio supernova ever detected. With a 6 cm maximum flux density of 
1.8 mJy, SN 1988Z is ^20% more luminous than the unusually powerful radio supernova SN 1986 J in 
NGC 891 and only ~3 times less radio luminous at 6 cm peak than the extraordinary SN 1998bw, the 
presumed counterpart to GRB 980425. Our analysis and model fitting of the radio light curves for SN 
1988Z indicate that it can be well-described by a model involving the supernova blastwave interacting 
with a high-density circumstellar cocoon, which consists almost entirely of clumps or filaments. SN 1988Z 
is unusual, however, in that around age 1750 days the flux density begins to decline much more rapidly 
than expected from the model fit to the early data, without a change in the absorption parameters. We 
interpret this steepening of the radio flux density decline rate as due to a change in the number density of 
the clumps in the circumstellar material (CSM) without a change in the average properties of a clump. If 
one assumes that the blastwave is traveling through the CSM at ~ 2, 000 times faster than the CSM was 
established (20,000 km s _1 vs. 10 km s -1 ), then this steepening of the emission decline rate represents a 
change in the presupernova stellar wind properties ~ 10,000 yrs before explosion, a characteristic time 
scale also seen in other radio supernovae. Further analysis of the radio light curves for SN 1988Z implies 
that the SN progenitor star likely had a ZAMS mass of ~ 20-30 M©. We propose that SNe, such as 
SN 1986J, SN 1988Z, and SN 1998bw, with very massive star progenitors and associated massive wind 
(M > 10 -4 Mq yr _1 ) have very highly-clumped, wind-established CSM and unusually high blastwave 
velocities (> 20,000 km s _1 ). 

Subject headings: supernovae, individual (SN 1988Z), galaxies, individual (MCG +03-28-022, Zw 
095-049), radio continuum, stellar evolution: massive stars 



1. INTRODUCTION 

Supernova (SN) 1988Z was independently discovered in 
MCG +03-28-022 (Zw 095-049) near m B ~ 16.4 by both 
G. Candeo on 1988 December 12 (Cappellaro, Turatto, & 
Candeo 1988) and C. Pollas on 1988 December 14 (Pollas 
1988). The SN was classified as a Type II, based on the 
detection of hydrogen emission lines in the optical spec- 
tra by Cowley & Hartwich (Heathcote et al. 1988), which 
also indicated that the SN was quite distant, with red- 
shift z ~ 0.022. Filippenko (1989) confirmed the Type 



II identification, with possible resemblance to SN 1987F. 
The spectra revealed at early times that SN 1988Z was 
peculiar, with a remarkably blue color, a lack of absorp- 
tion lines and P-Cygni profiles, very narrow (<100 km s _1 
FWHM) [O III] circumstellar emission lines, and a steadily 
growing, relatively narrow (~2000 km s _1 FWHM) com- 
ponent to the H I and He I lines (Stathakis & Sadler 1991; 
Filippenko 1991a,b). Stathakis & Sadler (1991), and Tu- 
ratto et al. (1993) for later times, analyzed in detail the 
spectroscopic and photometric observations of SN 1988Z. 
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sified as a "hypernova," approaching the 2-5 x 10 52 ergs 
estimated to have been released in SN 1998bw, the pos- 
sible counterpart of GRB 980425 (Iwamoto et al. 1998), 
and perhaps indicative of the collapse of the stellar pro- 
genitor core into a black hole. (It is interesting to note here 
that two SNe Iln, 1997cy and 1999E, may also have been 
associated with 7-ray bursts; see Pastorello et al. 2002 
and references therein). Aretxaga et al. (1999) suggest an 
ejecta mass of ~ 15 M Q , and, therefore, a very high- mass 
progenitor. 

Obviously, SN 1988Z is an extremely interesting and, 
in many ways, unusual object. Fortunately, due to its 
intrinsic brightness it has been relatively well-studied in 
many wavelength bands. Here we report radio observa- 
tions of SN 1988Z at multiple wavelengths, including new 
observations which add another six years of monitoring, 
more than doubling the coverage reported by Van Dyk et 
al. (1993a). Wc further interpret its radio emission using a 
clumpy wind model. We conclude that SN 1988Z, like SNe 
1986J and, possibly, other SNe Iln, as well as the unusual 
SN 1998bw, arise from the explosions of very massive stars 
surrounded by highly filamentary CSM. 

2. RADIO OBSERVATIONS 

New radio observations of SN 1988Z have been made 
with the VLA at 20 cm (1.425 GHz), 6 cm (4.860 GHz), 
3.6 cm (8.440 GHz) and 2.0 cm (14.940 GHz) from 1993 
May 4 through 2001 January 25 and are presented in Ta- 
ble 1 and Figure 1. The previously published results from 
Van Dyk et al. (1993a), along with a few additional and 
reanalyzed flux density values, are also included in Table 
1 and Figure 1 for ease of reference. 

Note that the estimated explosion date of 1988 January 
23 obtained by Van Dyk et al. (1993a) from their model 
fit was very poorly constrained by their fitting procedure. 
Here, we adopt the explosion date of 1988 December 1 es- 
timated from optical data (Stathakis & Sadler 1991). We 
have correspondingly modified the SN age at each of the 
epochs in Table 1 and Figure 1 relative to those reported 
in Van Dyk et al. (1993a). This change in assumed ex- 
plosion date, although large, is not critical for the model 
fitting and does not significantly affect the quality of the 
fit or the conclusions made by Van Dyk et al. (1993a). 

The techniques of observation, editing, calibration, and 
error estimation are described in previous publications on 
the radio emission from SNe (see, e.g., Weiler et al. 1986, 
1990, 1991). The "primary" calibrator was 3C 286, which 
is assumed to be constant in time with flux densities of 
14.45, 7.42, 5.20, and 3.45 Jy at 20, 6, 3.6, and 2 cm, re- 
spectively. The "secondary" calibrators 3 were J1051+213, 
used through 1991 January 17, and J1125+261, for all 
epochs after that, with defined positions of RA(J2000) 
= 10 h 51 m 48!789, Dcc(J2000) = +21°19'52'.'31 and 
RA(J2000) = ll h 25 m 53!712, Dec(J2000) = +26°10T9'.'98, 
respectively. After flux density calibration by 3C 286, 
they served as the actual gain and phase calibrators for 
SN 1988Z. As expected for secondary calibrators, the flux 
densities of J1051+213 and J1125+261 have been varying 

2 The VLA telescope of the National Radio Astronomy Observatory is operated by Associated Universities, Inc. under a cooperative agreement 
with the National Science Foundation. 

3 Secondary calibrators are chosen to be compact and unresolved by the longest VLA baselines. While compact and serving as good phase 
references, such objects are usually variable, so that their flux density must be recalibrated from the primary calibrators for each observing 
session. 



Schlegel (1990) proposed that SN 1988Z, along with other 
similar SNe, constituted a new class of SNe, the Type II- 
"narrow," or SNe Iln. 

The properties of the optical spectra and light curves 
indicated strong SN shock-circumstellar shell interaction 
(Filippenko 1991a,b; Turatto et al. 1993) which, to- 
gether with its resemblance optically to SN 1986J (Fil- 
ippenko 1991a, b; see also Rupen et al. 1987, Leibundgut 
et al. 1991), strongly suggested that SN 1988Z should be 
a very luminous radio emitter, as SN 1986J is a strong ra- 
dio source (see Weiler, Panagia, & Sramek 1990). Sramek 
et al. (1990) reported the radio detection at 6 cm wave- 
length of SN 1988Z with the Very Large Array (VLA) 2 on 
1989 December 21. The supernova is located at RA(J2000) 
= 10 h 51 m 5O!138, Dec(J2000) = +16°00'00'.'16, with an 
uncertainty of ± 0'.'2 in each coordinate, which is coinci- 
dent, to within the uncertainties, with its optical position 
(Kirshncr, Leibundgut, & Smith 1989). This initial an- 
nouncement of the detection of radio emission was followed 
by a more thorough study and analysis of the first three 
years of multifrcquency radio observations by Van Dyk et 
al. (1993a). They compared SN 1988Z in detail with SN 
1986J and concluded that the two SNe are very similar in 
their radio properties. They suggested that the progeni- 
tor to SN 1988Z was a massive [20 < M(M@) < 30] star 
which underwent a high mass- loss phase (M > 10 4 M Q 
yr~ 4 ) before explosion (see also Stathakis & Sadler 1991). 

Luminous radio emission has also been an indicator of 
X-ray emission from SNe, with synchrotron radio emission 
being produced as the forward SN shock interacts with 
the CSM and X-rays being emitted from the correspond- 
ing reverse shock region interacting with the SN ejecta 
(Chevalier & Fransson 1994). For instance, SN 1986 J 
was detected as a luminous X-ray source by Bregman & 
Pildis (1992) and Houck et al. (1998). Fabian & Terlevich 
(1996) reported the detection of X-rays from SN 1988Z 
with ROSAT. SN 1988Z is a very luminous X-ray emitter 
with a bolometric X-ray luminosity of ~ 10 41 erg s _1 . 

Chugai & Danziger (1994) offered two models to explain 
the unusual characteristics of SN 1988Z: 1) shock inter- 
action with a two component wind consisting of a ten- 
uous, homogeneous medium with embedded higher den- 
sity clumps, or 2) shock interaction with a similar tenu- 
ous, homogeneous medium and a higher-density, equato- 
rial, wind-established, disk-like component, favoring the 
former over the latter. However, they unexpectedly con- 
clude that the SN ejecta is of low mass (M < 1 M ), and 
that SN 1988Z may have originated from a relatively low- 
mass star of Mzams ~ 8-10 M Q , in sharp contrast to the 
high-mass progenitor suggested by Van Dyk et al. (1993a) 
and Stathakis & Sadler (1991). 

Aretxaga et al. (1999) collect the observations from X- 
ray to radio and attempt to estimate the integrated elec- 
tromagnetic energy radiated by SN 1988Z in the first 8.5 
years after discovery. They obtain a value of > 2 x 10 51 
erg, perhaps as high as 10 52 erg, which they consider is 
sufficiently high to suggest that SN 1988Z could be clas- 
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over the years, as can be seen in Table 2 and Figure 2. 

The flux density measurement errors for SN 1988Z are 
a combination of the rms map error, which measures the 
contribution of small unresolved fluctuations in the back- 
ground emission and random map fluctuations due to re- 
ceiver noise, and a basic fractional error e, included to 
account for the normal inaccuracy of VLA flux density 
calibration (see, e.g., Wcilcr ct al. 1986) and possible de- 
viations of the primary calibrator from an absolute flux 
density scale. The final errors (cr/) given for the measure- 
ments of SN 1988Z are taken as 



a) = (eSo) 2 + ol 



(1) 



where So is the measured flux density, cto is the map rms 
for each observation, and e = 0.05 for 20, 6, and 3.6 cm, 
and e = 0.075 for 2 cm. 

3. PARAMETERIZED RADIO LIGHT CURVES 

Following Weiler et al. (1986, 1990) and Montes, Weiler 
& Panagia (1997), we adopt a parameterized model (N.B.: 
The notation is extended and rationalized here from pre- 
vious publications. However, the "old" notation of r, t', 
and t", which has been used previously, is noted, where 
appropriate, for continuity): 



''"external 



^internal 



(2) 



with 



r external — T CSM hon 



where 



^CSM homogeneous 



t-tg Y 

1 day J 

2.1 



and 



rcsM clumps =r' = ^ 3 (^) (^) 



6' 



(1 



^clumps ) T( _ 



term describes the attenuation 



CSM c i umps 

produced by an inhomogeneous medium ( "clumpy absorp- 
tion" ; see Natta & Panagia 1984 for a more detailed dis- 
cussion of attenuation in inhomogeneous media). Both 
terms have a radial dependence which, under a constant 
mass-loss rate, constant velocity wind assumption, is r~ 2 
(see, e.g., Van Dyk et al. 1994 for an example where, for 
SN 1993J, the radial dependence of the CSM density is 
found to be flatter than r~ 2 ). The values of 5 and 5' de- 
termine the actual radial density profile if a constant shock 
velocity is assumed. The Tdistant term describes the atten- 
uation produced by a homogeneous medium which com- 
pletely covers the source, but is so far from the SN progen- 
itor that it is not affected by the expanding SN blastwave 
and is constant in time. All absorbing media are assumed 
to be purely thermal, singly ionized gas, which absorbs via 
free-free (f-f) transitions with frequency dependence v~ 21 
in the radio. The parameters 5 and 5' describe the time de- 
pendence of the optical depths for the local homogeneous 
and clumpy or filamentary media, respectively. 

The f-f optical depth outside the emitting region is pro- 
portional to the integral of the square of the CSM den- 
sity over the radius. Since in the simple model (Chevalier 
1982a, b) the CSM density decreases as r~ 2 , the external 
optical depth will be proportional to r~ 3 , and since the ra- 
dius increases as a power of time, r oc t m (with m < 1; i.e., 
m = 1 for undecelerated blastwave expansion), it follows 
that the deceleration parameter, m, is 



m 



-6/3. 



(7) 



(3) 

(4) 
(5) 

,5 GHz/ Vlday/ ' ^ 

with Ki, K 2, Ks, and K4 corresponding, formally, to the 
flux density (K\), homogeneous {K 2l K4), and clumpy 
or filamentary (K 3 ) absorption at 5 GHz one day af- 
ter the explosion date, to- The terms rcsM homo eneouB 
and TcsM c i umpB describe the attenuation of local, homo- 
geneous CSM and clumpy CSM that are near enough to 
the SN progenitor that they are altered by the rapidly 
expanding SN blastwave. T"cSM homogeneouB is produced by 
an ionized medium that homogeneously covers the emit- 
ting source ("homogeneous external absorption"), and the 

4 Note that, for simplicity, we use an internal absorber attenuation of the form 



The Chevalier model relates (3 and S to the energy spec- 
trum of the relativistic particles 7 (7 = 2a — 1) by 
5 = a — f3 — 3, so that, for cases where K 2 = and 5 
is, therefore, indeterminate, we can use 



m = —(a — ft — 3)/3. 



(8) 



Since it is physically realistic and may be needed in 
some RSNe where radio observations have been obtained 
at early times and high frequencies, we have also in- 
cluded in equation (2) the possibility for an internal ab- 
sorption term 4 . This internal absorption (Ti nt0 rnai) term 
may consist of two parts - synchrotron self-absorption 
(SSA; TintomaissA ) j an d mixed, thermal f-f absorption/non- 
thermal emission (Ti ntCT nai ff ) , so that 



^"internal — 'HnternalssA ^mternalff , 



where 



TmtcrnalssA 



U GHz/ ^ 



t-t 
1 day 



and 



^"intcrnalff 



GH2 



-2.1 



1 day / 



(9) 



(10) 



(11) 



which is appropriate for a plane-parallel 



T CSM intornal 

geometry, instead of the more complicated expression (e.g., Osterbrock 1974) valid for the spherical case. The assumption does not affect the 
quality of our analysis because, to within 5% accuracy, the optical depth obtained with the spherical case formula is simply three-fourths of 
that obtained with the plane-parallel slab formula. 
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with K§ corresponding, formally, to the internal, non- 
thermal (v a ~ 2 - 5 ) SSA and K 6 corresponding, formally, to 
the internal thermal {v~ 2A ) f-f absorption mixed with non- 
thermal emission at 5 GHz one day after the explosion 
date, to- The parameters 8" and 5"' describe the time de- 
pendence of the optical depths for the SSA and f-f internal 
absorption components, respectively. 

A cartoon of the expected structure of an SN and its sur- 
rounding media is presented by Weiler, Panagia, & Montes 
(2001; their Figure 1). The radio emission is expected to 
arise near the blastwave (Chevalier & Fransson 1994). 

The success of the basic parameterization and modeling 
has been shown by the good correspondence between the 
model fits and the data for all subtypes of RSNe, e.g., the 
Type lb SN 1983N (Sramek, Panagia, & Weiler 1984), the 
Type Ic SN 1990B (Van Dyk et al. 1993b), the Type II 
SNe 1979C (Weiler et al. 1991, 1992a; Montes et al. 2000) 
and 1980K (Weiler et al. 1992b, Montes et al. 1998), and 
SN 1988Z (Van Dyk et al. 1993a). (Note that, after day 
~ 4000, the evolution of the radio emission from both SNe 
1979C and 1980K deviates from the expected model evo- 
lution, and SN 1979C shows a sinusoidal modulation in its 
flux density prior to day ^4000.) A more detailed discus- 
sion of SN radio observations and of modeling results is 
given in Weiler et al. (2002). 

Thus, the radio emission from SNe appears to be rel- 
atively well understood in terms of blastwave interac- 
tion with a structured CSM, as described by the Cheva- 
lier (1982a, b) model and its modifications by Weiler et 
al. (1986, 1990) and Montes et al. (1997). For example, 
the fact that the homogeneous external absorption expo- 
nent 5 ~ — 3, or somewhat less, for most RSNe is evidence 
that the absorbing medium is generally a wind with den- 
sity p oc r~ 2 , as expected from a massive stellar progenitor 
which explodes during the red supergiant (RSG) phase. 

Additionally, in their study of the radio emission from 
SN 1986J, Weiler et al. (1990) found that the simple 
Chevalier (1982a, b) model could not describe the rela- 
tively slow turn-on. They therefore included terms de- 
scribed mathematically by TcsM c i umps in equations (2) and 
(6). This modification greatly improved the quality of the 
fit and was interpreted by Weiler et al. (1990) to represent 
the possible presence of filaments or clumps in the CSM. 
Such a dumpiness in the wind material was again required 
for modeling the radio data from SN 1993 J (Van Dyk et 
al. 1994) and SN 1988Z (Van Dyk et al. 1993a). Since 
that time, evidence for filamentation in the envelopes of 
SNe has also been found from optical and UV observations 
(e.g., Filippenko, Matheson, & Barth 1994: Spyromilio 
1994). 

Through use of this modeling a number of physical prop- 
erties of SNe can be determined from the radio observa- 
tions. One of these is the mass-loss rate from the SN pro- 
genitor prior to explosion. From the Chevalier (1982a, b) 
model, the turn-on of the radio emission for RSNe provides 
a measure of the presupernova mass-loss rate to wind ve- 
locity ratio (M/ui w ind). Weiler et al. (1986) derived this 
ratio for the case of pure, external absorption by a homo- 
geneous medium. However, we now recognize that several 
possible origins for absorption exist and generalize equa- 
tion (16) of Weiler et al. (1986) to 



^ (Me yr " } = 3.0 x 10- < > m-" x 



(w w i n d/10 km s *) 

\ I k 



10 4 km s _1 



45 days 



1.5 / , \ 1.5m 

k 



T 



10 4 K 



(12) 



Since the appearance of optical lines for measuring SN 
ejecta velocities is often delayed relative to the time of 
explosion, we normally adopt t\ = 45 days. Because our 
observations generally have shown that 0.8 < m < 1.0, 



and from equation (12) M oc t i 



1.5(l-m) 



, the dependence of 



the calculated mass-loss rate on the date ti of the initial 
ejecta velocity measurement is weak, M oc tf . Thus, 
we generally adopt the best optical or VLBI velocity mea- 
surements available, without worrying about the deviation 
of their exact measurement epoch from the assumed 45 
days after explosion. For SN 1988Z we follow Turatto et 
al. (1993) and assume Vi = 20, 000 km s _1 from the highest 
velocity measured for the broad components of the emis- 
sion lines. One has to keep in mind that the actual shock 
velocity for SN 1988Z may be even higher, because it is 
very hard to derive accurate velocities from fitting very 
broad lines where the extreme wings may be lost in the 
continuum. 

We also normally adopt values of T = 20, 000 K, 
Wwind = 10 km s _1 (which is appropriate for a RSG wind), 
t = {tg cm peak — to) days from our best fits to the radio 
data, and m from equation (7) or (8), as appropriate. The 
optical depth absorption term, < r^j 5 >, however, is not 
as simple as that used by Weiler et al. (1986). 

Weiler et al. (2001) were able to identify at least three 
possible absorption regimes: 1) absorption by a homoge- 
neous external medium, 2) absorption by a clumpy or fil- 
amentary external medium with a statistically large num- 
ber of clumps, and 3) absorption by a clumpy or filamen- 
tary medium with a statistically small number of clumps. 
Each of the three cases requires a different formulation for 
< Tgg 5 >. From our consideration of the radio light curves 
for SN 1988Z, we conclude that the Case 3 of Weiler et 
al. (2001, 2002) is most appropriate (see §4). 

Case 3 occurs when the clump number density is small, 
the probability that the line-of-sight from a given clump 
intersects another clump is low, and both emission and 
absorption effectively occur within each clump. Such a 
situation will yield a range of optical depths, from zero 
for clumps on the far side of the blastwave-CSM interac- 
tion region, to a maximum corresponding to the optical 
depth through a clump for clumps on the near side of the 
blastwave-CSM interaction region. We expect an atten- 
uation of the form, (1 — e ~ TCSMclumps ) T csM c i umps ' as gi yen 
in equation (2), but now TcsM clumpB represents the optical 
depth along a clump diameter. Since the clumps occupy 
only a small fraction of the volume, they have a volume 
filling factor <j> <C 1. Additionally, the probability that the 
line-of-sight from a given clump intersects another clump 
is low, so that a relation between the size of a clump, the 
number density of clumps, and the radial coordinate can 
be written as 

r, Trr 2 R w N < 1, (13) 
where r\ is the volume number density of clumps, r is the 
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radius of a clump, R is the distance from the center of the 
SN to the blastwave-CSM interaction region, and N is the 
average number of clumps along the line-of-sight, with N 
appreciably lower than unity by definition. Finally, it is 
easy to verify that there is a relation between the volume 
filling factor (j>, r, R and N, of the form 



4 r 
3 R 



N. 



(14) 



We can then express the effective optical depth < > 

as 



< r eff >- 



V2 



0.5 

r CSM c 



0.5 j\ r 0.5 



(15) 



where, for initial estimates, we shall take N ~ 0.5 and a 
constant ratio ri? _1 ~ 0.33, so that, from equation (14), 
<p ~ 0.22. Before we can fully estimate the mass-loss rate, 
M, for the progenitor star of SN 1988Z before explosion, 
we must first fit our parameterized radio light curve model 
to the radio dataset for the SN. 

4. RADIO LIGHT CURVE DESCRIPTION 

Examination of Figure 1 and Table 1 shows that the 
multi-frequency radio data are best described by two time 
intervals: "early" data, which extends roughly from ex- 
plosion through day 1479, and "late" data, which roughly 
extends from day 2129 through the end of the dataset. Fig- 
ure 1 shows clear steepening of the light curves sometime 
between these two measurement epochs, but the actual 
"break" date is somewhat arbitrary, due to uncertainties 
in the flux density measurements for this relatively faint 
source and the likely smoothness of any transition region. 
We have therefore chosen to describe the flux density evo- 
lution separately for these two time intervals, with the pe- 
riod from day <~ 1500 to day <~ 2000 as a transition. With 
the explosion date assumed to be 1998 December 1 from 
optical estimates (Stathakis & Sadler 1991), and exten- 
sive trial fitting showing no evidence for distant, homoge- 
neous thermal absorption (K4 = 0), for SSA (K§ = 0), for 
mixed free- free absorption/nonthermal emission (Kq — 0), 
or for a local homogeneous thermal absorption component 
(K 2 = 0), the model fitting process needs to determine 
only five parameters (K\, a, 0, K$, and 5'). Applying our 
fitting procedures separately to the early and late periods, 
with the data points between day 1500 and day 2000 in- 
cluded in both fits to provide a smooth transition, we find 
that both the spectral index, a, and the clumpy absorp- 
tion parameters, K3 and 5' , are the same in the two time 
intervals, to within the fitting errors. However, the emis- 
sion decay rate parameter, 0, steepens significantly, from 
<~ — 1.2 for the early period to ~ —2.7 for the late period. 
The results of this model fitting are listed in Table 3. Note 
that K\ is a flux density scaling factor, so the fact that it 
has very different values for the early and the late periods 
is not physically significant. 

For a purely clumpy CSM (K 2 = 0), the sharp 
steepening of around day 1750, while K3 and 5' re- 
main unchanged, implies that: 1) the number density of 
clumps per unit volume, 77, starts decreasing more rapidly 
with radius by approximately i? -1 - 5 /™ (i.e., ryd a y>i750 = 
-1.5/n 

%ay<1750( flT " 



), with the average characteris- 



tics of the individual clumps remaining unchanged, and 2) 
most of the absorption occurs within the emitting clumps 
themselves. In other words, the spatial distribution is so 
sparse that the average number of clumps along the line- 
of-sight is less than one (N < 1). This second condition is 
consistent with Case 3 from Weiler et al. (2001, 2002) and 
is the basis for our selection of Case 3 for SN 1988Z. It is 
perhaps significant that Weiler et al. (2001, 2002) found 
that Case 3 also applies to the radio light curves for the 
unusual SN 1998bw. 

Considering other wavelengths, we note that the Ha 
light curve from Aretxaga et al. (1999; see their Figure 2) 
deviates significantly from their models after day — 1000. 
To investigate this further, we have plotted in Figure 3 
the Ha light curve after day 300, the interval for which 
we have radio data. Examination of Figure 3 shows that 
the Ha light curve is consistent with a slope steepening 
from 1 = -1.84 to j3 2 = -2.94 after day -1250. The 
Ha emission is proportional to p 2 , while for "typical" SNe 
the radio emission is proportional to p 1A (i.e., if j3 is the 
Ha light curve decline rate, 0\ = —1.84 should be the 
same as — 2/1.4 x /3 ra dio = 1-42 x /3 oar i y = —1.73, where 
Nearly = —1-22 from Table 3). Thus, very good agreement 
is found between the behavior of the decline for both the 
early optical Ha emission-line and early radio light curves, 
to within the uncertainties. Also, the change in the slope 
(A/3) at both Ha and in the radio should be the same, be- 
cause the break is due to a change in the number of clumps, 
not in the properties of the clumps and, therefore, should 
be wavelength independent: A/? H a = (-2.94) - (-1.84) = 
-1.10; A/? radio = (-2.73) - (-1.22) = -1.51, which is also 
in fair agreement, considering the uncertainties. Thus, it 
appears that the steepening in the decline for the nonther- 
mal radio light curves is consistent with the steepening of 
the thermal Ha light curve. 

4.1. The Mass-Loss Rate for the SN Progenitor 

We can now use equations (12) and (15) to estimate from 
the radio absorption a mass-loss rate for the SN 1988Z pro- 
genitor star. With the assumptions for the blastwave and 
CSM properties discussed in §3 and the results for the 
best-fit parameters listed in Table 3, we obtain an esti- 
mated presupernova mass- loss rate, M = 1.2 x 10~ 4 M Q 
yr" 1 . 

This high mass-loss rate for SN 1988Z is only appropri- 
ate for the last — 10,000 years before explosion (see §5). 
At earlier epochs of the progenitor's evolution the mass- 
loss rate was considerably lower, as indicated by the much 
more steeply declining (3 with unchanged K% and 6' after 
day -1750. 

5. DISCUSSION 

Figure 1, Table 3, and §4 show that the radio light curves 
for SN 1988Z can be described by standard RSN models 
(Weiler et al. 1986, 1990, 2001; Montes et al. 1997), with 
only one fitting parameter varying with time — the in- 
dex of the time evolution of the radio emission, (3. At an 
age of ~1750 days (3 clearly steepened from — —1.2 to 
— —2.7. This change in radio emission characteristic, 0, 
was not accompanied by a corresponding change in the 
absorption characteristics, K 3 and 5'. Thus, the change in 
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the radio light curves for SN 1988Z cannot be described as 
only due to a change in the average CSM density, which 
so effectively described the changes in radio evolution for 
SN 1980K (Montes et al. 1998) and SN 1979C (Montes et 
al. 2000). 

For SN 1988Z the observations are best interpreted as 
a decrease in the number density of clumps in the wind- 
established CSM after day ~ 1750, but with all clumps 
having essentially the same average internal characteristics 
as at earlier times. Since the shock velocity for SN 1988Z 
is assumed to be ~ 20, 000 km s _1 and the pre-supernova 
wind velocity for an RSG is typically ~ 10 km s _1 , an in- 
terval of 1,750 days after explosion implies that this change 
took place in the presupernova stellar wind ~ 9, 600 years 
before explosion. With an estimated mass-loss rate in the 
interval prior to explosion of <~ 1.2 x 10~ 4 M yr -1 , it 
follows that during this last 10,000 years, the progenitor 
star shed ~ 1.1 M Q in a massive "superwind." The abrupt 
steepening of f3 around day 1750 implies that the mass- loss 
rate was appreciably lower before that time, decreasing by 
about a factor of four by day 4,438. Even such a decreas- 
ing mass-loss rate, however, still accounts for an additional 
~ 0.8 Mq lost in the stellar wind by our last measurement 
on day 4,438 (~ 25, 000 years before explosion, using the 
same wind and blastwave velocity assumptions). Since a 
relatively high mass-loss rate must have been maintained 
over a much longer time, such as the ~ 10 5 years appro- 
priate for a massive star's evolution, an additional several 
Mq must have been lost earlier so that the presupernova 
mass loss for SN 1988Z over the entire RSG progenitor 
lifetime is expected to have been > 2 Mq, and perhaps 
much greater. 

As discussed in §1, Chugai & Danziger (1994) studied 
SN 1988Z in detail and developed a model from which they 
concluded that the ejecta mass is small (M e j < 1 Mq), 
and that SN 1988Z, therefore, must have originated from 
a relatively low-mass star, Mzams ~ 8-10 Mq, in sharp 
contrast with the high-mass (20-40 Mq) progenitor pro- 
posed by Van Dyk et al. (1993a) and Stathakis & Sadler 
(1991). However, it should be noted that Chugai & 
Danziger derived a very high mass-loss rate for the SN 
1988Z progenitor, 7 x 10~ 4 M yr -1 , yielding what they 
call their "wind density parameter" (our M/w w i n d) of 
5 x 10 16 g cm -1 . We believe that our mass-loss rate de- 
rived from the direct measurement of f-f absorption in 
the CSM, 1.2 x 10~ 4 M & yr _1 (wind density parameter, 
M/w w i nc i = 7.3 x 10 15 g cm -1 ), is likely more appropri- 
ate. Since the critical factors in the modeling by Chugai & 
Danziger are the product of the ejecta mass and wind den- 
sity parameter related to the energy of the explosion, for a 
given explosion energy our 6-7 times smaller wind density 
parameter implies, for their model, a factor of 6-7 times 
greater ejecta mass, M e j ~ 5-10 Mq, rather than their es- 
timate of M e j < 1 Mq. Even our value for the ejecta mass 
may still be an underestimate because, by the Chugai & 
Danziger model, it represents only the mass that was ob- 
served to interact with the pre-supernova wind and may 
not involve the entire ejecta mass. With an ejecta mass of 
as much as ~ 10 Mq and a high radio luminosity, which 
places SN 1988Z among the brightest radio supernovae, 
we conclude that the stellar progenitor of SN 1988Z must 
have been a very massive star, perhaps as high as ~ 20 - 



30 M Q , as previously suggested by Van Dyk et al. (1993a). 

It is worth noting that all of the well-studied cases 
of ultra-bright RSNe (say, having >10 times the radio 
luminosity of SN 1979C at 6 cm peak), such as SNe 
1988Z, 1986J, and 1998bw, provide evidence for a highly- 
clumped CSM with presupernova wind mass-loss rates up 
to <~ 1O~ 4 M yr" 1 and unusually high shock velocities 
> 20,000 km s _1 . From optical and ultraviolet spec- 
troscopy Fransson et al. (2002) also find for the SN Iln 
1995N evidence for high shock velocities, a "superwind" 
phase for the progenitor characterized by a high mass-loss 
rate prior to explosion, and a clumpy CSM. Such high- 
rate mass loss, clumping of the CSM, and consequent high 
blastwave velocities may be the signature of particularly 
massive star progenitors (see also the discussion in Frans- 
son et al. 2002). It is yet to be determined how the progen- 
itors of SNe Iln, which must retain at least some of their 
hydrogen envelopes until the end of their lifetimes, and 
the progenitors of possibly extreme Type Ib/c SNe, such 
as SN 1998bw, which must have lost their entire hydrogen 
envelopes before explosion, are related, if at all. 

6. CONCLUSIONS 

The radio emission from SN 1988Z followed an evolu- 
tion well described by standard models up to an age of 
~ 1750 days (~ 4.8 yrs), after which its behavior evolved 
into a much more rapid decline in radio emission without 
a corresponding change in radio absorption parameters. 
This is the first time we have witnessed this in any well- 
studied RSN, and it implies that highly radio-luminous 
RSNe, such as SN 1988Z and, possibly, SN 1986J, may 
follow a different evolutionary path in their presupernova 
mass-loss than that for more "normal" RSNe, such as SNe 
1979C and 1980K. However, it is interesting to note that, 
if one assumes a blastwave velocity ~ 2, 000 times faster 
than the RSG wind-established CSM (~ 20, 000 km s" 1 
vs. ~ 10 km s~ 4 ), the timcscale for such a variation in SN 
1988Z of - 10 4 yr, is similar to that seen for SNe 1979C 
and 1980K (see, e.g., for SN 1979C, Montes et al. 1998; for 
SN 1980K, Montes et al. 2000), and also for SN 1998bw 
(Weiler et al. 2001). 

In addition, further analysis of the results from Chugai 
& Danziger (1994), which implied an extremely high mass- 
loss rate of 7 x 10~ 4 M Q yr -1 and a relatively low-mass pro- 
genitor with Mzams ~ 8-10 M for SN 1988Z, indicates 
that their estimates are probably unrealistic. Using their 
modeling, our lower mass-loss rate, 1.2 x 10 Mq yr -1 , 
implies a significantly higher ejecta mass and, therefore, a 
higher ZAMS mass of ~ 20-30 Mq for the progenitor star. 

We have noted that the Ha data from Aretxaga et 
al. (1999), which are also indicative of the SN shock-CSM 
interaction (Chevalier & Fransson 1994), also begin to de- 
viate significantly from their model by day >1250. Al- 
though a detailed comparison is beyond the scope of this 
paper, both the slope before the break and the magnitude 
of the break are roughly consistent with the radio results. 

From our analysis we propose that SNe, such as SNe 
1988Z, 1986J, and 1998bw (possibly the counterpart to 
GRB 980425), with possible very massive star progeni- 
tors (Mzams > 20 Mq) and associated massive winds 
(M > 10~ 4 Mq yr -1 ), have very highly-clumped, wind- 
established CSM and unusually high blastwave velocities 
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(> 20,000 km s- 1 ). 
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Table 1 



Flux Density Measurements for SN 1988Z a 


Obs. Date 


Age 


VLA 


& (20 cm) 


o (6 cm) 


S (3.6 cm) 


S (2.0 cm) 


(UT) 


(days) 


Config. 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


1988 Dec 1 


= 












1989 Dec 21 


385 


D 




0.67 ±0.07 






1990 Feb 12 


438 


DnA 




0.77 ±0.09 






1990 May 29 


544 


A 




1.21 ±0.09 






1990 Jul 12 


588 


BnA 




1.26 ±0.09 


2.10 ±0.14 




1990 Sep 7 


645 


B 


<0.39 


1.38 ±0.09 


1.72 ± 0.14 


1.15 ±0.23 


1990 Dec 14 


743 


C 


<0.48 


1.78 ± 0.10 


2.09 ±0.13 


1.43 ±0.19 


1991 Jun 11 


922 


A 


0.47 ±0.07 


1.90 ±0.11 


1.68 ± 0.10 


1.31 ±0.25 


1991 Sep 17 


1020 


A 




1.85 ± 0.10 


1.48 ±0.09 


0.68 ±0.16 


1992 Jan 26 


1151 


CnB 


1.04 ±0.13 


1.67 ±0.10 


1.50 ±0.09 


1.17 ±0.18 


1992 Oct 13 


1412 


A 


0.91 ±0.08 


1.57 ±0.09 


1.22 ± 0.10 


0.91 ±0.22 


1992 Dec 19 


1479 


A 


1.52 ±0.09 


1.72 ± 0.11 


1.04 ±0.08 


0.78 ±0.19 


1993 May 4 


1615 


B 


1.16 ±0.10 


1.21 ±0.09 


0.98 ± 0.06 


0.57 ±0.15 


1993 Aug 23 


1726 


BnA 


1.05 ±0.15 


1.28 ± 0.09 






1993 Dec 5 


1830 


D 






0.87 ±0.07 




1994 Feb 18 


1905 


DnA 


1.09 ±0.16 


0.98 ± 0.17 


0.86 ±0.14 




1994 May 4 


1980 


BnA 


1.42 ±0.16 


0.97 ±0.13 


0.52 ± 0.08 




1994 Sep 30 


2129 


CnB 


1.06 ±0.15 


0.91 ±0.05 


0.56 ± 0.08 




1995 May 19 


2360 


D 




0.62 ± 0.05 


0.42 ± 0.04 


0.26 ±0.09 


1995 Oct 3 


2497 


BnA 


0.78 ±0.07 


0.40 ± 0.04 


0.23 ±0.04 


<0.38 


1996 Oct 5 


2865 


DnA 


0.72 ±0.15 


0.36 ± 0.09 


0.32 ± 0.06 


<0.36 


1997 Jan 23 


2975 


BnA 


0.68 ±0.11 


0.30 ±0.08 






1998 Feb 10 


3358 


DnA 






< 0.15 




1998 Feb 13 


3361 


DnA 


0.35 ±0.06 


0.30 ±0.05 






1998 Oct 19 


3609 


B 


0.39 ± 0.07 


0.20 ± 0.03 






1999 Jun 13 


3846 


DnA 


0.19 ±0.06 


<0.20 






1999 Oct 1 


3956 


DnA 


0.22 ±0.06 


0.17 ±0.04 






2001 Oct 19 


4340 


A 


0.19 ±0.03 


<0.18 


<0.18 




2001 Jan 25 


4438 


DnA 


0.17 ±0.06 


<0.26 


<0.27 





a All upper limits are 3<r. 
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Table 2 

Flux Density Measurements for Secondary 
Calibrators 



Obs. Date S (20 cm) S (6 cm) S (3.6 cm) S (2.0 cm) 
(UT) (Jy) ' (Jy) (Jy) (Jy) 



J1051+213 



1989 Dec 21 




0.921 






1990 Feb 12 




0.885 






1990 May 29 




0.901 






1990 Jul 12 




0.907 


0.894 




1990 Sep 7 


1.092 


1.002 


0.942 


1.123 


1990 Dec 14 


1.118 


0.948 


1.069 


0.967 


1991 Jun 11 


0.962 


0.997 


1.162 


1.203 


1991 Sep 17 




1.091 


1.358 


1.336 



J1125+261 



1992 Jan 26 


0.992 


1.069 


1.001 


0.806 


1992 Oct 13 


0.948 


1.102 


0.954 


0.839 


1992 Dec 19 


0.937 


1.068 


0.926 


0.729 


1993 May 4 


0.967 


1.008 


0.871 


0.703 


1993 Aug 23 


0.892 


0.986 






1993 Dec 5 






0.887 




1994 Feb 18 


0.910 


1.049 


0.934 




1994 May 4 


0.920 


1.056 


0.891 




1994 Sep 30 


0.943 


1.082 


0.984 




1995 May 19 




1.088 


0.936 


0.749 


1995 Oct 3 


0.880 


1.098 


0.966 


0.781 


1996 Oct 5 


0.899 


1.059 


0.854 


0.657 


1997 Jan 23 


0.893 


1.060 






1998 Feb 10 






0.918 




1998 Feb 13 


0.961 


0.911 






1998 Oct 19 


0.989 


0.967 






1999 Jun 13 


0.957 


0.927 






1999 Oct 1 


0.937 


0.988 






2000 Oct 19 


0.924 


0.887 


0.709 




2001 Jan 25 


0.933 


1.070 


0.936 
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Table 3 

Fitting Parameters for SNe 1988Z 



Parameter 




SN 1988Z a 






Early 
day < 2000 


Global 


Late 
day > 1500 


a 

K 3 
5' 

h h 

L 6 cm peak (erg s 1 Hz 1 ) 

t = (*6 cm peak - to) (days) 

M {Mq yr-^ f 


1.2 x 10 4 
-1.22 


-0.72 

3.19 x 10 8 
-2.87 
eee 1988 Dec 01 
2.3 x 10 28 

911 
1.2 x 10~ 4 


9.1 x 10 8 
-2.73 



a Becausc of a change in the evolution of the radio emission from SN 
1988Z during the interval between day ^1500 and day ~2000, the data are 
split into two overlapping intervals of "early," from day 385 through day 
1980, and "late," from day 1615 through day 4438. 

b Adoptcd from Stathakis & Sadler (1991). 

c Assuming w wind = 10 km s~\ Vi = fbiastwave = 20,000 km s _1 , T = 
20, 000 K, ti = 45 days, N = 0.5, rR- 1 = 0.33, and cf) = 0.22 and calculating 
to = 0.83 and < t°J> >= 0.16. 
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Fig. 1— Radio "light curves" for SN 1988Z in MCG +03-28-022. The four wavelengths, 2 cm (14.9 GHz, open circles, solid line), 3.6 cm 
(8.44 GHz, crosses, dashed line), 6 cm (4.87 GHz, open squares, dot-dash line), and 20 cm (1.48 GHz, open triangles, dotted line), are shown 
together with their best-fit model light curves. The SN age is in days since the adopted explosion date of 1988 December 1 (Stathakis & 
Sadler 1991). Because the decline index f3 of the radio emission steepened between, roughly, day ~1500 and day ^2000, an iterative fitting 
procedure was used (see §4), and a break in the model radio light curves can be seen near day 1750. 



Fig. 2. — Flux density measurements of the secondary calibrators J1051+213 and J1125+261 at 2 cm (14.9 GHz, open circles), 3.6 cm (8.44 
GHz, crosses), 6 cm (4.87 GHz, open squares), and 20 cm (1.48 GHz, open triangles). For clarity the data points at each frequency are offset 
from one another by adding 1.0 Jy to the 2 cm values, 0.5 Jy to the 3.6 cm values, and —0.5 Jy to the 20 cm values. The epoch at which the 
secondary calibrator was changed from J1051+213 to J1125+261 is marked. 
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Fig. 3. — Evolution of the Ha line intensity, from Aretxaga et al. (1999), during the same time range when radio observations were obtained, 
~300 to ~3000 days after explosion. Note a likely steepening of the decline rate, f3, for the Ha intensity after ~ 1250 days, with j3\ = —1.84 
before the break and /32 = —2.94 thereafter. 



This figure "fl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0208190vl 



This figure "f2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0208190vl 



